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Dark energy dynamics in the recent universe is influenced by its evolution through the long, 
matter dominated expansion history. A particular dynamical property, the flow variable, remains 
constant in several classes of scalar field models as long as matter dominates; the dark energy is 
only free to diverge in behavior at recent times. This gives natural initial conditions for Monte 
Carlo studies of dark energy dynamics. We propose a parametrization for the later evolution that 
covers a wide range of possible behaviors, is tractable in making predictions, and can be constrained 
by observations. We compare the approach to directly parametrizing the potential, which does not 
take into account the maturity of the dark energy dynamics. 



I. INTRODUCTION 

Very little is known about the dynamics of dark en- 
ergy, other than that at redshifts z > 1 the expansion 
history gives way to a matter dominated era. However, 
this fact is powerful in narrowing the possible behav- 
iors of dark energy because the matter-dominated Hubble 
friction plays a major role in the dark energy dynamics. 
That is, the dark energy field <f> evolves for a long time 
in a particular environment - by the time dark energy 
comes to contribute significantly to the expansion it is 
an old, or "mature" , field. 

The influence of matter domination on the dark energy 
dynamics is very different than the slow-roll conditions 
for an inflationary field in the early universe. Dark en- 
ergy does not satisfy slow-roll conditions on the potential 
V for most of its evolution, unless it is highly fine tuned 
Instead, the matter domination can create a defi- 
nite relation between the deviation of the dark energy 
density fi^ from 0, the deviation of its equation of state 
w from —1, and the characteristic scale of the potential, 
(l/V)dV/d<p (which is sometimes called the first slow-roll 
parameter, but here does not need to be small). A partic- 
ular combination of them is called the flow parameter by 
and is essentially constant during matter domination. 

In this paper we examine the dark energy evolution af- 
ter matter domination wanes, by parametrizing the devi- 
ation of the flow parameter from its constant value dur- 
ing matter domination. Since the constancy holds up 
until quite late redshifts, z « 2, the dark energy has 
a restricted ability to exhibit diverse dynamics by the 
present. This allows the three parameter flow form we 
present to reasonably approximate a wide range of the 
possible behaviors. 

Any form containing few parameters will not in general 
be able to describe every possible behavior so interesting 
questions include whether it captures key physics, reason- 
ably describes the major classes of behavior, and can be 
constrained by observations. In particular, the dynamics 
should include the thawing and freezing classes [3] , which 
respectively depart from and approach cosmological con- 
stant behavior within specific regions of the w-w' phase 



space. But it should also be flexible enough to permit 
models outside these regions and have the physics and 
data determine which are viable. 

In Sjn]we discuss the implications of a long matter dom- 
inated era in guiding the dark energy dynamics and the 
role of the flow parameter. We present our approach to 
evolution as dark energy increases in importance in ijllll 
and study the resulting behaviors. Constraints on the 
w-w' and flow parameters phase spaces from future ob- 
servations are investigated in iflVI We compare the flow 
approach to parametrizing the potential in SjVj discussing 
the different weights these impose on the dynamics and 
the thawing vs. freezing classes. 

II. MATURE DYNAMICS 

During the early universe accelerated expansion of 
inflation, the inflaton scalar field dominates the total 
energy density and so the field evolution is essentially 
wholly determined by its potential. Once inflation ends, 
the cosmic expansion is dominated by radiation and then 
matter. The expansion rate - the Hubble parameter H - 
is rapid and the Hubble drag influences the evolution of 
any light scalar field that may eventually cause a more 
recent universe acceleration. 

The evolution of this scalar field is determined by the 
Klein-Gordon equation 

$ + 3H{p b ,pt)<j> + V <l> = Q, (1) 

where we have explicitly indicated the dependence of 
the Hubble parameter on the background energy density 
Ph and the scalar field density p^. The energy density 
P0 = V(4>) + (l/2)</> 2 , the sum of the potential and ki- 
netic energies, so we cannot determine the field evolution 
without knowing V((f>) and the initial conditions <j>i and 
<j)i. Formally, then, we cannot say anything completely 
general about the dark energy dynamics. 

However, an initial kinetic energy much larger than the 
potential, K 3> V, will rapidly redshift away, <p ~ a~ 3 as 
can be seen by neglecting the last term in Eq. ([T]). So we 
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expect the kinetic energy of the field does not dominate, 
but rather the background fluid (radiation or matter) 
does. See 0-H for early papers on cosmological scalar 
fields, or quintessence. The distance A(f> the field rolls can 
be estimated by approximating (j) rj Ac/)/ At and using 
K = (1/2) p^, (1 + w), where the equation of state w — 
[(K/V)-1))/[(K/V) + 1]. Then 



A0~ yj(l + w) Pl p (At) 2 



(2) 



If the Hubble friction is large, then the characteristic 
timescale of the field evolution may be the Hubble time 
H , so the square root contains p^/H 2 ~ f2^ which is 
small in the matter (or radiation) dominated era. So the 
era of background domination has significant effect on 
the dark energy dynamics. 

If the potential dominates the kinetic energy then 1+w 
is very small as well and so both factors lead to A<fi <C 1 
during the background domination (although by today, 
with dark energy domination, the condition A(f> <C 1 gen- 
erally does not hold; see Sec.fV]). This class is known as 
thawing fields. In other cases 1 + w does not start small 
and the field can roll further. 

If the potential has particular forms, the field can pos- 
sess an attractor trajectory (see, e.g., d, 04I3l)- In 
this case the equation of state is determined by the back- 
ground component and the parameters of the potential, 
independent of initial conditions. The potential and ki- 
netic energies are locked in a set ratio such that the equa- 
tion of state is constant during background domination; 
the fields are then known as trackers. 

In any case, the dark energy density is not exactly zero 
and the equation of state is not exactly —1. It is interest- 
ing to consider how these deviations are related to each 
other in the background dominated era. Ref. Q showed 
(also see [ll|) that one could define a combination of 
these deviations, together with the characteristic scale of 
the potential, that held constant during the long matter 
dominated expansion. Ref. [2j defined 



F = 
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(3) 



where 1+w is the equation of state deviation, the dark 
energy density in units of the critical density, and A = 
— {\/V)dV/d(f> gives a characteristic field scale. Despite 
all these quantities varying as the universe expands and 
the field evolves, the flow combination F remains almost 
constant. 

Written in terms of the dark energy dynamics, 



w 



-3(1 



1 



1 



(4) 



where w' — dw/d\na. Thus the long matter dominated 
era, which imposes F — constant until quite late times, 
z « 2, (except for specifically fine tuned fields) effectively 
imposes constraints on the dark energy dynamics. That 
is, the dark energy is an old field that lived for many 



e-folds of expansion in a matter dominated environment. 
Choosing an arbitrary dynamics amounts to either ignor- 
ing this physical influence or embracing a fine tuning of 
the field's initial conditions. 

A direct parametrization of the potential (see, e.g., 

EMI), 

and Monte Carlo variation of its parameters 
a la the interesting dynamics work of [HI, [13], will not 
in general take into account the maturity. While one 
could parametrize the equation of state in such a way 
as to incorporate the appropriate condition (see, e.g., 
[EH EL EH), we choose to study the flow parameter since 
it possesses the interesting properties of being nearly con- 
stant during matter domination and incorporating O^, 
the natural quantity for measuring the deviation from 
matter domination as it begins to break. 



III. EVOLVING DYNAMICS 



During matter domination, fields which slowly relax 
from being frozen by the high Hubble friction - thawing 
fields - possess F = 4/27 Q- Equivalently, they thaw 
from the cosmological constant state along the trajectory 
w' = 3(1 +w) according to Eq. (U). Tracker fields possess 
F = 1/3 or w' = 0. Fields with other behaviors are 
certainly possible, but these behaviors will be sensitive 
to the specific initial conditions. However in this section 
we will briefly consider other high redshift values for F. 

We are most interested in how the dark energy evolves 
as matter domination wanes, and as eventually dark en- 
ergy comes to dominate by the present. That is, we want 
to extend the asymptotic, matter dominated behavior 
to greater dynamical freedom in the present. Calculat- 
ing the first order (in the dark energy density fraction) 
corrections to the dynamics Q gives a flow parameter 
behavior 



F w F + (3 a 1 



(5) 



where Fq is the matter dominated constant value, /3 is 
related to the dark energy density fi^, and B = 3 (or 
—3 Woo where w(z 3> 1) = Woo in the thawing (or track- 
ing) case. 

However, since today the ratio of dark energy density 
to matter density is nearly 3, we do not expect the first 
order solution to be valid. Numerical solutions, for ex- 
ample, show that F can grow by factors of order two to 
three. Therefore, although we will keep the form ([5]) we 
take /3 and B to be free parameters representing effec- 
tive fits over the expansion history between high redshift 
matter domination and today. 

The form |JSJ| has the advantage that, upon substitution 
into Eq. Q, one can analytically solve for the dynamics 
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w(a). The result is 

x-i 
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where Fi — F(ai), Wi — w(ai), and is some initial 
scale factor, e.g. in the matter dominated regime. The 
four parameters are therefore F 0l /3, B, and w i7 and we 
choose a% = 0.1 (the exact value of aj is unimportant as 
long as it is well within the matter dominated era). 

As the flow F increases as the dark energy density 
becomes more important, the field evolves further in the 
direction of the skating curve, w' — — 3(1— w 2 ), where the 
potential slope is unimportant and field moves - skates 
[Til [l9| - with a rapidly redshifting kinetic energy. For 
thawing fields, this means that the dynamics can actually 
turn around and head back toward the frozen, cosmolog- 
ical constant state. Recall that thawing fields start with 
Fq = 4/27, so their maximum, turnaround value of w 
occurs when F has risen to 1/3, at which point w' = 0. 
For freezing fields, the increase in F simply moves them 
faster toward the cosmological constant. 

Figure [TH5] show the effects of each parameter on the 
dynamics, for the thawing and freezing cases. The ampli- 
tude parameter /3 controls the degree of evolution from 
the matter dominated state - increasing (3 has the same 
effect as increasing F. The transition parameter B gov- 
erns the rapidity of evolution - larger B means that the 
transition from the matter dominated behavior happens 
nearer to the present. 

The initial condition parameters Fq and Wi have a 
significant influence, as seen in Fig. [3] Increasing Fq 
gives the dynamics a significant push in its evolution. 
If Fq 7^ 4/27 then the field somehow ignores the high 
Hubble friction during matter domination. Basically the 
field must have been "super-accelerating" , i.e. <f> ^> Hep, if 
Fo < 4/27, while it must have been fine tuned to be coast- 
ing 4> <C H(j) (not slow-roll since the matter, not the dark 
energy, dominates the expansion) if 4/27 < F < 1/3 
(cf. [l| ) . The value F = 1/3 corresponds to the dynam- 
ics determined by an attractor solution for the field, but 
F > 1/3 indicates "super-deceleration", — » dV/dcj), 
equivalent to fine tuning for a "super-flat" potential, yet 
somehow the field still rolls. These properties motivate 
why we will later consider only the values Fq — 4/27 and 
1/3. 

Similarly, increasing Wi for thawing fields gives a strong 
boost to the evolution, pushing the dynamics far from 
w = — 1. (For tracker freezer fields, Wi is determined by 
the attractor solution and then = 0.) For thawing 
fields, we expect 1 + Wi to be of order fi^j, the dark 
energy density at the initial redshift considered in the 
matter dominated era (see, for example, [2]). (Taking 
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FIG. 1: [Top panel] The deviation amplitude of the flow pa- 
rameter between the matter dominated era and today is given 
by p. Increasing /3 moves the dynamics away from small field 
deviations and toward skating, with rapid approach toward 
a frozen state. [Bottom panel] For a fixed deviation ampli- 
tude, the transition rapidity is determined by B. Increasing B 
postpones the transition from the matter dominated flow pa- 
rameter, causing a more rapid evolution close to the present. 
Both panels are for the thawing case, with Fo = 4/27, and 
I + Wi — 10 -4 ; curves extend to a = 2 with thick portions of 
curves for a < 1. 



di = 0.1, typically ~ - 1 + 10 4 .) For freezing fields Wi 
can be of order unity, but again tracking evolution starts 
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FIG. 2: As the bottom panel of Fig. [1] but for the freezing 
case, with Fo = 1/3, and Wi = —0.5. Variation of /3, as in 
Fig.[T]top panel, tends to move the evolution along the curves; 
for example for B = 3 we show by x's the a = 1 location for 
/3 = 1(T 5/2 , 10~ 2 , 1(T 3/2 , 10~\ 10" 1/2 , from right to left. 



with w' — 0. Another aspect of non-standard values of 
Fo or Wi is the presence of excessive early dark energy 
density. For example, a thawing field with 1 + Wi — 
1CP 2 has fi^j ps 0.06, which would cause conflict with 
growth of large scale structure and the cosmic microwave 
background observations [2(| . A thawing field with Fo = 
1.5/27 would break matter domination completely. 

Thus, we effectively have a three parameter descrip- 
tion, given by the set {(3,B,Wi}, for each of two cases: 
F) = 4/27, 1/3. In general even a three parameter de- 
scription of dark energy dynamics cannot be fully con- 
strained even with next generation data [Il|,[22j]. Indeed, 
determining the transition time or speed for the dark en- 
ergy equation of state is extremely difficult [2l|, [23|, [24| 
so we do not expect B to be constrained. However, this 
form is useful because it gives considerable freedom for 
dynamics while at the same time taking into account the 
oldness of dark energy, i.e. the long influence of matter 
domination on the evolution. 

The ansatz also provides for three classes of dynami- 
cal behavior: the usual thawing and freezing classes, and 
"re-freezing" models where the field thaws from the Hub- 
ble friction-induced torpor, evolves to a maximum w, but 
then can enter the freezing region and eventually evolve 
to a potential minimum and cosmological constant state. 
This last class will be especially interesting as it bounds 
the area of the w-w' plane that lies within a certain ob- 
servational distance of A, i.e. models that agree with the 
ACDM expansion history to within a certain precision 
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FIG. 3: [Top panel] The flow parameter during matter dom- 
ination is generally determined to be Fo = 4/27 or 1/3, but 
other values are possible with fine tuning. Non-standard 
values can cause unusual, strong dynamical behavior. For 
Fo < 6.5/27 we show curves taking Wi = —0.99; larger values 
of Fo give unobservably small loops. For Fo > 1/3 we take 
Wi — —0.5. [Bottom panel] The dark energy equation of state 
at some fixed redshift in the matter dominated era can give a 
strong thrust to the evolution for thawing models if it is not 
near the thawing value u>i = — 1 expected from Hubble drag. 



(see the next section). One can imagine a fourth class 
which starts as trackers but then turns away from the 
cosmological constant, i.e. somehow moves away from the 
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minimum, but there is little physical motivation for this 
and we do not consider it further. 



IV. CONSTRAINING DYNAMICS 

In the previous section we saw that rich dynamics is 
still available to the dark energy field, even after taking 
into account its oldness - the effect of the long era of 
matter domination to age and mellow the arbitrariness 
of the early evolution. In particular, we emphasized the 
naturalness of the age influence and did not impose any 
restrictions on behavior such as assuming small field ex- 
cursions all the way until late times. Here we consider 
how observational data will be able to constrain the later 
time behavior within the flow parameter formulation and 
narrow in on regions of the w-w' plane. 



A. Scanning Phase Space 

In fact, we will find clear, bounded regions of the phase 
space and maximum allowed deviations in 1 + w, as- 
suming some future distance data. To scan over all the 
possible dynamics within the flow ansatz we begin by 
considering the parameter ranges. As discussed in the 
previous section, there are good reasons from both the- 
ory and observations to consider initially thawing and 
tracker/freezer classes, with F = 4/27 and 1/3 respec- 
tively. We emphasize that these values only give the ini- 
tial conditions deep in the matter dominated era and are 
not assuming later behavior. 

For 1 + Wi , we indicated that for thawing models this 
was often around 10~ 4 , while for freezing models it could 
lie between a small value and 1, though it is more dif- 
ficult to achieve w rj — 1 today if Wi is not sufficiently 
negative. However, to keep open all possibilities, we con- 
sider 1 + Wi ranging between 10 -5 and 1, with a log prior 
for the thawing class. For the amplitude /?, most models 
generally considered have values between fewxlCP 2 and 
fewxlO -1 . Again to be flexible, we take a log prior be- 
tween 10 -5 ' 2 and 1. The evolution rapidity B can be ex- 
amined from the first order correction (in the small dark 
energy density) to matter domination; for thawing mod- 
els B = 3 and for tracking models B — —3wi. Numerical 
solutions show that for many thawers B slowly declines 
toward the present, reaching ~ 1.5, and for trackers B 
slowly rises toward the present, reaching ~ 2. We take a 
generous linear prior on B between 1 and 5. 

The models obtained by scanning over the parameter 
space on a 50 x 50 x 50 grid in (3-B-Wi, for F — 4/27 and 
1/3, are then compared in the distance-redshift relation 
to the distances to redshifts z = 0.1, 0.2,. . . 1.7 taking a 
ACDM fiducial cosmology. The 250,000 generated mod- 
els are referred to as the prior sample, and those that pass 
observational cuts are referred to as the viable sample. 
The observational cut considered first is simply requiring 
the distances at each redshift to lie within 1% of fiducial 



(roughly next generation data precision), and the dark 
energy density at m = 0.1 must obey n^j < 0.03. We 
later consider a more sophisticated Monte Carlo likeli- 
hood calculation. 

Figure 0] shows the dynamics in terms of the evolution 
of the equation of state with scale factor, w(a), for the 
prior and viable models. While the observations certainly 
constrain the viable models to more limited values of the 
equation of state, the variety of dynamics allowed by the 
flow parametrization is still diverse. In particular, models 
with thawing (leaving w = —1), freezing (approaching 
— 1), and refreezing (nonmonotonic in w) behavior are 
all represented. 




FIG. 4: The flow ansatz allows a rich variety of evolutionary 
behavior, including nonmonotonicity, while incorporating the 
physical influence of the long matter dominated era. The top 
panel shows a random selection of 150 models (note many 
are very close to w = — 1) while the bottom panel shows 150 
random models that pass the observational cuts. (Note the 
different vertical scales.) Thawing, freezing, and nonmono- 
tonic evolution are all represented. 

Observational constraints translate to well defined re- 
gions in the w-w' phase space, with a maximum allowed 
distance deviation imposing a maximum allowed equa- 
tion of state deviation w max . Furthermore, the allowed 
region is compact, limiting w' as well. Figures [5][6] shows 
the regions that have less than 1% distance deviation 
from the ACDM fiducial, with the colors/shading repre- 
senting different values of Wi, as well as scanning over 
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j3 and B. Note that in Fig. [5] for small 1 + Wi most 
models look like regular thawing models, while as l + Wi 
increases, the only way the models can obey the observa- 
tional constraints are to loop around as refreezing models. 
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FIG. 5: Models deviating from ACDM distances by less than 
1% live in a well defined region of phase space, here shown for 
the initially thawing class. Colors/shading represent different 
values of Wi , with models possessing the smallest values of 1 + 
Wi restricted to the bright, red region closest toui = —1. The 
top panel shows models with Wi < —0.999 (as expected for 
most thawing models); the bottom panel shows models with 
Wi > —0.999. The sparseness of the bottom panel is due to the 
finite gridding: only a tiny percentage of the 2500 models for 
each Wi pass the distance criteria. The maximum equation of 
state deviation allowed at any redshift corresponds to w max = 
-0.845. 



B. Mirage Model Dynamics 

An interesting special case to consider is those mod- 
els that deviate from ACDM by exactly 1%, rather than 
those in the whole range 0-1%. These could represent, for 
example, a specific fiducial distinguishable from ACDM 
by next generation data. Figures [7][8] show the w(a) evo- 
lutionary behaviors. Of particular note is the waist near 
a = 0.8, where the equation of state has a dispersion of 
only Aw « 0.01. 

This sort of waist feature has been discussed before in 
the context of mirage models [25[ and pivot points, albeit 
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FIG. 6: As Fig. [5] but for the initially freezing class. The 
maximum initial equation of state deviation, at en — 0.1, 
corresponds to w max = —0.36, reaching Wo = —0.985 today. 

restricted to the parametrization w(a) = wo + w a (l — a). 
Such a characteristic highlights the importance of exper- 
iments capable of measuring the time variation w' of the 
equation of state; a survey only able to measure a coarse 
time average will ineluctably find (w) ~ 0.95 for all these 
diverse models. If the dynamical model allowed the equa- 
tion of state to become more negative than —1, then one 
could easily obtain a waist with ~0% distance deviation 
from ACDM and find (w) — —1, wrongly deducing a 
cosmological constant as the answer. This is the mirage 
aspect of the mirage models, and so uncovering the true 
physics behind dark energy requires an experiment with 
a long redshift baseline over which to measure the time 
variation of the equation of state. 



C. Dynamics Today 

We can also consider the properties of the dark energy 
today. For initially thawing models, the values of wq and 
w' can lie basically anywhere within the overall observa- 
tionally allowed region, e.g. the top panel of Fig. [5] (for the 
bottom panel, where Wi > —0.999, the values today lie 
at the end of the loops, within the top panel region). For 
initially freezing models, the distribution today is more 
interesting. For Wi near —1, the values today fill most 
of the appropriate part of the freezing region in Fig. [5] 
The region allowed by 1% distance observations is largest 
when Wi ~ —0.95. As Wi increases, however, the dynam- 
ics today hugs the lower boundary of the freezing region, 
w' = 3w(l+w), and the allowed region shrinks, vanishing 
for Wi > —0.36. Figure [9] illustrates these characteristics. 

To understand the elements of the dynamics, we show 
in Fig. [10] the role of the parameters /? and B on the 
present value of the equation of state and its time varia- 
tion, holding Wi fixed. As expected, /3 is more influential 
than B, with large values of /3 accelerating the evolution 
toward a present state near the cosmological constant. 

The original classification of the boundaries of the 
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FIG. 7: Initially thawing models that have maximum distance 
deviations from ACDM by some fixed amount, here 0.999- 
1%, have a variety of behaviors but exhibit a narrow waist at 
a « 0.8. In terms of a constant or averaged equation of state, 
these will all look like (w) ~ —0.95 but have distinct physics. 
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FIG. 8: As Fig. but for freezing models. Note the waist 
persists, though here it is at a w 0.7. 

thawing and freezing regions Q was made by studying 
numerical solutions for a range of field potentials and 
determining their endpoints, i.e. wo-w' values today. Al- 
though the long period of matter domination guides the 
flow approach dark energy models toward an initial be- 
havior either along the upper boundary of the thawing 
region, w' = 3(1 + w), or to the tracker condition of 
w' = and subsequent approach toward the cosmologi- 
cal constant, we can now address quantitatively whether 
their evolution leads to consistency with the thawing and 
freezing regions today. 

We clearly see from Fig. [5] that the initially freezing 
models lie through most of their evolution in the canon- 
ical freezing region 0.2w(l + w) < w' < 3w(l + w), and 
from Fig. [5] that in particular they all lie in the freezing 
region today. Recall though that thawing models in the 
flow approach can either remain thawing, loop around, 
or even refreeze. Out of all the viable (0-1% distance 
deviation) initially thawing models, 47% remain in the 
thawing region through the present, while 40% today lie 
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FIG. 9: The dynamics today, wo-w' Q , are shown with each 
dot representing the present endpoint of a trajectory in Fig. [6] 
for freezing models giving distance deviations of 0-1% from 
ACDM. As Wi moves away from —1, the region gradually 
grows, reaching a maximum extent at Wi ~ —0.95 (top panel), 
then shrinking and coalescing along the bottom boundary of 
the freezing region (delimited by the straight lines). The bot- 
tom panel shows the case of Wi = —0.8. 



in what would usually be the freezing region, and 13% 
today lie between the thawing and region regions near 
the coasting (4> = 0) line. If we restrict consideration 
to the most natural initially thawing models, those with 
1 + Wi < 5 x 10 -4 (recall most standard thawers have 
1 + Wi w 10 -4 ), then the percentages become 60%, 28%, 
12% respectively. Thus there is no shortage of thawing 
models today. 



D. Monte Carlo Constraints 

In the previous subsections we made use of a grid tes- 
sellation of parameter space to examine the variety of 
w-w' phase space behaviors that can be obtained within 
the flow formalism. Here we will probe constraints on the 
F(a) parameter space. In particular, we are interested in 
the ability to distinguish between model classes given ob- 
servational data, and secondarily to determine the values 
of the flow parameters. We focus on the capability of fu- 
ture observations to rule out the cosmological constant, 
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FIG. 10: For the present values of the equation of state pa- 
rameters wo and w' , the influence of flow parameters f3 and B 
are illustrated. The isocontours of /3 (thick lines) are shown 
for ,9 = 1, 0.11, 0.033, 0.01, 0.0032 from bottom to top for 
the top panel (/3 = 1, 0.5, 0.2, 0.11, 0.05, 0.02 from left to 
right for the bottom panel). Isocontours of B (thin lines) 
are shown for B = 1, 2, 3, 4, 5 from top to bottom. The 
top panel shows the initially thawing case, the bottom panel 
the initially freezing case, with the dashed lines showing the 
conventional boundaries for each region. 



for example through measuring an evolution through a 
nonzero /? or finding the high redshift (z > 10) equation 
of state value Wi different from — 1. 

To quantify this, we generate data using a fidu- 
cial model and perform a Markov Chain Monte Carlo 
(MCMC) likelihood analysis [26M30lj on the F(a) param- 
eters which, together with the initial value of the equa- 
tion of state parameter Wi, fully determine the evolution 
of the dark energy. The usual ACDM model corresponds 
to = —1 and the best fit is achieved with a set of 
parameters that select from the prior the smallest value 
of 1 + Wi (to start close to A), the largest of f3 (to make 
the field quickly freeze, i.e. approach A), and the smallest 
value of B (to start the freezing evolution early). Basi- 
cally the cosmological constant behavior is an asymptotic 
state for scalar field dynamics. Because of this, the like- 
lihood would concentrate around the prior bounds for 
whatever ranges we choose, making the predictions prior 
dependent. To avoid this, we instead choose a fiducial 



cosmology distinct from ACDM. This is of further inter- 
est since it explores the sensitivity of F(a) parameters in 
revealing this deviation, should next generation datasets 
narrow in on a region of parameter space off ACDM. 

Following the approach of the previous sections we 
study thawing and freezing models separately since these 
have well defined, high-redshift responses to the flow 
function F(a). Accordingly, we select Fq equal to 4/27 
and 1/3, respectively. As discussed in Sec. Mil fields with 
initial F values different from these either betray matter 
domination, dilute very fast, or are finely tuned coasters. 
Nonetheless we emphasize that Fq merely determines the 
high redshift initial conditions of the field and other pa- 
rameters govern its late time dynamics. The distinction 
between thawing and freezing origin motivates further 
a different choice of priors for 1 + Wi, since the differ- 
ent initial behavior of w(a) is one of the features distin- 
guishing between the two classes of models. For thawers, 
1 + Wi is related to the initial energy density fi^j and 
can be quite small, so we choose a logarithmic prior on 
1 + Wi for the thawing models; for trackers, 1 + Wi tends 
to be larger so we choose a uniform prior on 1 + Wi for 
the (tracker) freezer case. Similarly, since f3 tends to be 
larger in the freezing case we choose a uniform prior there 
but a logarithmic prior for the thawer case. The bounds 
on the priors, however, are the same as for the previous 
grid scanning: 1 + Wi € [1CT 5 ,1], (3 € [1CT 5/2 ,1], and 

b e [1,5]. 

The fiducial cosmology for both classes of models is 
chosen to be a maximum distance deviation of ~ 1% 
away from ACDM, to test the constraint capability of 
next generation observations. We assume data of 1% 
accuracy in distance at redshifts z = 0.1, 0.2,. . . 1.7. This 
allows us to explore the degree to which the data can 
distinguish between the fiducial (thawing or freezing), its 
opposite class (freezing or thawing), and the cosmological 
constant. When considering a thawing fiducial we adopt 
F = 4/27, 1 + w t = 1Q- 4 , (3 = 0.04, and B = 2. For a 
freezing fiducial we take F = 1/3, Wi = —0.75, j3 = 0.6, 
and B = 1.5. 

The MCMC sampling of parameter space is done with 
the Metropolis-Hastings algorithm, using a Gaussian pro- 
posal distribution and fours chains per case. Computa- 
tion of probability distributions and likelihood contours 
is done using a modified version of the GetDist pack- 
age, provided with the CosmoMC software [28|. We de- 
fine Wi = w(ai) with cii = 0.1 and define the present by 



n 



0,0 



0.72. We integrate the coupled equations 



w' = -3(l-w) 2 



1 - 



1 



n'. 



d' 



-3u>r^ (i - n^) 



(1 - ri^o) a" 



(9) 
(10) 

(11) 



where d(a) is the conformal distance and a prime denotes 
d/d In a. The likelihood is based on the comparison of 
d(a) for a F(a) model to d(a) for the fiducial. 
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The probability distributions for the thawing and 
freezing fiducial cases are shown in Fig. [11] (left and 
right panels respectively). We find W{ to be the pa- 
rameter best constrained, with clear distinct regions of 
excluded parameter space. In particular, the likelihood 
pulls away from the lower bound of the prior, preferring 
that 1 + Wi > 10~ 5 and thus demonstrating a distinc- 
tion from the cosmological constant. The flow deviation 
amplitude ft is less well constrained, but for the thaw- 
ing fiducial case the likelihood has some preference for 
smaller values, in contrast to the freezing class. For the 
freezing fiducial case, the contour prefers larger values of 
ft, disfavoring a fit by thawing models. As mentioned 
previously, the evolution rapidity B is essentially uncon- 
strained. We emphasize that the flow approach is not 
intended as a new parametrization to fit specific mod- 
els, but rather a useful tool for distinguishing classes of 
physics, in particular with respect to the high redshift 
behavior. 

For both thawing and freezing fiducial cases the 2- 
dimensional joint probability distributions for (wi, ft) and 
(wi,B) are each almost decorrelatcd. The exception is 
for large values of ft (large overall flow F(a)) where w(a) 
need not start as negative since the evolution is pushed 
strongly toward ACDM. For (wi,B) the correlation is 
even weaker, and Wi is insensitive to B for most of its 
range except for very small values of B (corresponding 
to an earlier flow deviation, evolving towards ACDM), 
where u>i can be further away from the cosmological con- 
stant behavior. These results show that, while as ex- 
pected a three parameter description of the equation of 
state cannot be tightly constrained, the key discriminat- 
ing parameter of Wi can be constrained. This allows 
distinction (for this 1% deviating fiducial at least) be- 
tween each of the main classes of dark energy: a static 
cosmological constant, an initially thawing field, and an 
initially freezing field. 



V. COMPARING APPROACHES 

Throughout most of this article we propounded the use 
of the flow quantity F as a natural and slowly-evolving 
function, which could be relied upon to be approximately 
constant over a large number of e-folds. Unlike the case 
of parametrization of the scalar field potential, for which 
we need to assume some knowledge of its form and hence 
restriction to a particular case, the flow function F(a) 
possesses natural characteristics that allow us some in- 
sight on the dynamics of the field, on the sole assumption 
of matter domination, without having to assume some 
explicit form for V{<f>). 

It is interesting to explore this argument further by 
showing that the opposite approach ; a Taylor expansion 
of the potential as is often used (cf. (MlHl), 

can lead to 

results which don't fully or faithfully describe the avail- 
able and viable dynamics. Indeed its consistency as a de- 
scription of the potential can break down, as the quanti- 



ties in which it is being expanded are not generally small. 
Such is the case of the (not necessarily slowly rolling) 
slow-roll parameters e(4>) and ??(</>), defined in terms of 
first and second derivatives of the potential, or the field 
displacement A<j>. 

The usual parametrization of the potential based on 
the slow-roll assumption takes 

v(<p)/v(<p = Q)~i + v 1 4> + v 2 <p 2 + ... (12) 

where <po (taken to be zero without loss of generality) 
is the field value at some initial scale factor at, Vq is 
the initial value of the potential energy density, and the 
coefficients of the expansion are expressed in terms of the 
usual slow-roll parameters [3ll - l36j with V\ = —^/Wirei, 
and V2 — 47rrjj. 

For the validity of Eq. (|T2"f as a Taylor expansion, one 
needs either to rely on a sufficiently flat potential, by 
having small slow-roll parameters initially, or to take the 
expansion to be valid over a narrow region of the field tra- 
jectory, such that for the range of potential we are inter- 
ested in (i.e. its evolution up to the present) A(f> is small. 
In particular, we would want A<f> <C (167re) _1 ^ 2 and 
Acj) <C (e/lnr) 2 }) 1 ^ 2 to have convergence, i.e. the linear 
being smaller than the constant term, and the quadratic 
term smaller than the linear term, respectively. 

However, remember there is no guarantee the potential 
is sufficiently flat, today or at high redshift. This in turn 
gives no grounds for taking small bounds on the slow-roll 
parameters, in Monte Carlo simulations or other analy- 
ses. Indeed the detailed analysis of [uj finds values of 

~ 1 and r\i ~ 5, or greater, to be within the 68% confi- 
dence limits contour at z% = 3 (for current distance data 
taken from supernovae, baryon acoustic oscillation, cos- 
mic microwave background, and Hubble constant mea- 
surements; note they use the potential as a quadratic or 
cubic form, not a Taylor expansion per se) . If instead one 
tries to ensure validity of a (truly) slow-roll expansion by 
artificially imposing small (or equivalently starting the 
field at rest), one obtains strongly restricted dynamics as 
seen in the clear analysis and Fig. 9 of [13], which shows 
much narrower w{a) evolutions than in the absence of 
this imposition. 

To highlight the influence of matter domination on the 
slow-roll parameters in the potential model, we tested 
starting the initial conditions at different redshifts and 
following the evolution of e and r\. If instead of imposing 
(y L = 1, rji = 5 at z,i = 3, we start with these values at 
higher redshift, Zi = 9, then by z — 3 the first slow-roll 
quantity has dropped to e ~ 0.01 (77 stays fairly con- 
stant), clearly demonstrating the governing influence of 
matter domination. 

This reduction is insensitive to the initial redshift - 
as long as it is at least one e-fold before the evaluation 
point at z = 3. Thus, if one wants to carry out Monte 
Carlo simulations, either one can use a wide range for the 
parameters - but the initial conditions should be set at 
least a few e-folds into matter domination, or for a more 



10 





FIG. 11: 1- and 2-dimensional marginalized probability distributions and contours of MCMC reconstruction of the flow pa- 
rameters are illustrated for next generation distance data. The left panel shows results for the thawing case, with Fo = 4/27, 
which used log priors on 1 + Wi and /3; the right panel shows the distributions for the tracker/freezing case, with Fo = 1/3, 
which used uniform priors to take into account the larger fiducial values. The diagonal panels correspond to the marginalized 
ID probability distributions and the off diagonal to the marginalized 2D joint distributions. Thick curves in the off diagonal 
plots show the 68% and 95% likelihood contours and the color grading denotes the marginalized probability density from (red) 
highest to (blue) lowest. The MCMC constraints clearly recognize small 1 + to, for thawing, distinct from the large values 
typical of freezing models and from the zero value of the cosmological constant. For freezing, again Wi is distinct from the 
cosmological constant and from typical thawing values. The parameter /3 shows tendencies toward the expected preference for 
small values for thawing and larger values for freezing, but the constraints are not tight. Most importantly, the thawing and 
freezing classes can be distinguished from each other and from the cosmological constant. 



recent start one should use a restricted range of initial 
conditions as given by the matter dominated physics. 

In addition, in order to apply the Taylor expanded po- 
tential formalism one needs to rely on a small field excur- 
sion during the relevant dark energy evolution. In terms 
of the slow-roll parameters, for e ~ 1, r) w 5 the con- 
vergence conditions mentioned above require Acj) <C 0.14 
and Acj> <C 0.11. With the aim of testing this assumption 
we compute the field displacement Acj) incurred for the 
flow models we considered previously, between the initial 
redshift and today, conservatively selecting those which 
most resemble ACDM, since A(j> is smallest for these. We 
will find that we cannot guarantee the necessary level of 
smallness to validate the potential approximation. 

From the expression for the scalar field kinetic energy 
(see above Eq. [5]), the effective field excursion Acj> be- 
tween the initial redshift, Zi = 9 and today is 

A0 = / T V 3 ^ (1 + w) ■ (13) 

We compute Acj) for the flow models from Sec. IIVI that 
pass the 1% distance criteria. For thawing models the 
main contribution to Acj> is during the last e-fold or so, 



and accordingly we find a relatively small field displace- 
ment for thawers, of Acf> w 0.06 for Wi = — 1 + 10~ 5 and 
up to Acj> w 0.4 for the top of the range Wi = —0.99. 

On the other hand for the freezing behavior there is 
the possibility of a more substantial contribution to the 
field excursion throughout the evolution since the dark 
energy density is larger at early times and the equation of 
state finds its way towards A from larger Wi. For freezers 
we obtain 0.014 < Acj> < 0.53 between the smallest to 
the largest Wi, —0.99 < w t < —0.36. 

Both for thawers and freezers the leading influence on 
Acj) among the flow parameters is Wi . For thawers this is 
straightforward since w(a) increases constantly with u>i 
- apart from those models which cross F — 1/3 and turn 
around, receding back toward ACDM: for those Acj> will 
be smaller than in the purely thawing case. 

The freezing dark energy case is different: large Wi will 
cause the field to travel longer initially, before it slows 
down towards w ~ — 1 in order to meet the distance 
criteria. For these fields particularly, it is crucial to ini- 
tialize the evolution equations at high enough redshift, 
where the field is evolving the most and determining the 
conditions for its later dynamics. Starting the evolution 
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at z = 3, only 1.4 e-folds before the present, is prob- 
lematic and can change the results. For example, using 
zi = 3 Ref. [13] finds an average of A</» = 0.09 ± 0.03, 
out of their posterior distribution obtained from current 
data, 74% of the models being freezer models. Our re- 
sults using Zi = 9 for the flow model give significantly 
different results; even for models that closely resemble A 
we find values as large as A<fi = 0.34 between z = 3 and 
z = 0, for freezer models. Interestingly, in the absence 
of the slow-roll assumption either in the form of small 
slope or small field displacement, one finds the allowed 
field dynamics to be both richer and simultaneously able 
to remain close to ACDM. 

This argument becomes all the more important if we 
remember we don't have any guarantee today of distance 
measurements agreeing with ACDM to 1%. Relaxing this 
would further release the behavior of freezer models early 
on, and can give large field excursions, for example A<f> = 
0.87 for an n = 1 SUGRA model [13], with a maximum 
3% distance deviation for z < 1.7. 

Thus the flow parameter approach appears to have a 
certain amount of natural motivation, from the long era 
of matter domination, and can offer a self-consistent ap- 
proach. It is also important to note that the appropri- 
ate initial conditions must be carefully considered for a 
Monte Carlo analysis of the field evolution. 

Because the relation w'(w) can be written either in 
terms of the flow F or the potential parameters, e.g. 
A = — (l/V)dV/dcf> and Cls, interchangeably through 
F = (1 + wj/lfl^X 2 ], the near constancy of F implies 
that one cannot arbitrarily choose A, fl^, and Wi in the 
Monte Carlo initial conditions during matter domination. 
Moreover, if any of 1 + w, or A 2 start too large (in- 
cluding taking a uniform prior out to a relatively large 
upper bound) , then generating realizations of the thaw- 
ing class of models will be strongly suppressed. As men- 
tioned previously, even a large initial value of e (~ A 2 ), 
if the evolution takes into account matter domination, 
becomes a much more restricted value by z — 3 (this ac- 
cords with Figure 9 of [l7| where once they restrict the 
initial value of e, a substantial fraction of randomly gen- 
erated models become otherwise scarce thawers). While 
the flow approach is not without flaw, it does present 
an interesting alternative, one that innately incorporates 
the role of the long matter dominated era and the oldness 
of dark energy. 

VI. CONCLUSIONS 

Dark energy does not exist in a vacuum. For dozens 
of e-folds of Hubble expansion it was dominated by mat- 
ter and radiation, and has risen above them for only a 
fraction of the last e-fold. This has a definite effect on 



the dynamical evolution of dark energy, aging or mellow- 
ing it. In particular, a wide range of models exhibit a 
nearly constant flow parameter, relating the deviation of 
the equation of state from —1, the deviation of the dark 
energy density from 0, and the slope of the potential. 
It is not until z > 2 that the flow parameter, and the 
equivalent relation between w' and w, is generally free to 
vary. 

This characteristic suggests the flow parameter is an 
interesting quantity to use for parametrizing the dark 
energy evolution, one that can intrinsically take into ac- 
count the oldness of dark energy. We have examined 
an evolutionary model for the deviations of the flow 
parameter from constancy as matter domination wanes 
and acceleration occurs. This gives rise to a variety 
of dynamical behaviors in the w'-w phase space, in- 
corporating thawing, freezing, and "refreezing" classes. 
Observational constraints on the dark energy equation 
of state, e.g. through distance-redshift measurements, 
bound specific regions in the w'-w plane. For exam- 
ple, scanning 250,000 models, we illustrate the region, 
and its component classes, of maximum deviation al- 
lowed away from w = — 1. We also identify an interest- 
ing "waist" phenomenon - independent of the standard 
Wo-Wa parametrization - where at a certain redshift the 
equation of state is tightly constrained. 

With an MCMC analysis one can quantify the bounds 
on the values of the flow parameters allowed by distance 
measurements. In particular this permits determination 
of the high redshift dark energy equation of state, a key 
parameter for the important physics question of distin- 
guishing between a cosmological constant, initially thaw- 
ing models, and initially freezing models. Finally, we em- 
phasize that taking into account the maturity of dark en- 
ergy is important for MCMC generation of models since 
random initial conditions, e.g. on the potential, do not 
appropriately sample the "diversity under constraint" of 
old dark energy. 
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